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Abstract 



Quantum gravity phenomenology opens up the possibility of probing Planck scale 
physics. Thus, by exploiting the generic properties that a semiclassical state of the com- 
pound system fermions plus gravity should have, an effective dynamics of spin-1/2 parti- 
cles is obtained within the framework of loop quantum gravity. Namely, at length scales 
much larger than Planck length ip ~ 10~ 33 cm and below the wave length of the fermion, 
the spin-1/2 dynamics in flat spacetime includes Planck scale corrections. In particular 
we obtain modified dispersion relations in vacuo for fermions. These corrections yield a 
time of arrival delay of the spin 1/2 particles with respect to a light signal and, in the 
case of neutrinos, a novel flavor oscillation. To detect these effects the corresponding par- 
ticles must be highly energetic and should travel long distances. Hence Neutrino Bursts 
accompanying Gamma Ray Bursts or ultra high energy cosmic rays could be considered. 
Remarkably, future neutrino telescopes may be capable to test such effects. This paper 
provides a detailed account of the calculations and elaborates on results previously re- 
ported in a Letter. These are further amended by introducing a real parameter T aimed 
at encoding our lack of knowledge of scaling properties of the gravitational degrees of 
freedom. 
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1 Introduction 



It is commonly accepted that quantum gravity should hold at scales near Planck length ip 



\/ ^Newton V° 3 ~ lCT 33 cm or, equivalently, Planck energy Ep := hc/£p ~ 10 19 GeV. Accord- 
ingly, neither astrophysical observations nor ground experiments were considered in the past as 
a means to directly reveal any quantum gravity effect but only to test indirect consequences. In 
recent years however this attitude has changed on the basis of potentially testable phenomena 
probing quantum gravity scenarios in which scales combine to lie not far from experimental 
resolution. Prominent among such phenomena we find in vacuo dispersion relations for gamma 
ray astrophysics |], ^, ]3|, f|, || 0, laser- interferometric limits on distance fluctuations |7], ||, neu- 



trino oscillations |4], |9|, threshold shift in certain high energy physics processes |K| [TTJ, |12, [L3|], 
CPT violation [14| and clock-comparison experiments in atomic physics |]T5[| . These are the 
prototypes of the emerging quantum gravity phenomenology ]TB| , |T7| , |T8| , . 

The present work is aimed at elaborating on effective corrections to propagation in vacuo for 
spin-1/2 particles in the framework of loop quantum gravity as reported in |3Q. This framework 
also has been used previously in studying light propagation ||, || . 

In essence, the specific structure of spacetime for a given quantum gravity scenario can be 
probed by matter propagating and interacting there: what could be considered as flat spacetime 
macroscopically, might produce microscopic imprints of its detailed structure on the interaction 
of particles. In particular, dispersion relations of propagating matter could exhibit corrections 
due to such effects. For particles with energy E <<c Ep and momentum p the following in 
vacuo modified dispersion relations were proposed @ : 

where Eqg ~ Ep and £ ~ 1. In general, such corrections might behave as ( > where 

T > 0, namely not necessarily an integer. This possibility has been considered for photons 
0, |21| and spin zero particles ||21|| . In this work we show it applies also to spin 1/2 particles. 
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According to Eq.([l[) the corresponding particle's speed yields a retardation time 

*»f#-f (2) 

with respect to a speed c signal after traveling a distance L. Interestingly, for Gamma Ray 
Bursts (GRB's) with E ~ 0.20MeV, L ~ 10 10 ly and setting E QG ~ E P , the naive estimation 
Eq. (HD gives At ~ 0.01ms, barely two orders of magnitude below the sensitivity St for current 
observations of GRB's |22|, ^3| . Moreover it is expected to improve this sensitivity in the 
foreseeable future |24|]. Using an expression analogous to Eq.(|2|) for the delay of two photons 
detected with an energy difference AE, the observational bound Eqg/C, > 4 x 10 16 GeV was 
established in Ref. [25]| by identifying events having AE = 1 TeV arriving to earth within the 
time resolution of the measurement At = 280 s, from the active galaxy Markarian 421. 

Now, accompanying GRB's there seems to be Neutrino Bursts (NB) in the range 10 5 — 10 10 
GeV according to the so-called fireball model |26], If detected they could provide an excellent 
means to test quantum gravity effects of the type given by Eq.(|2|) above. Experiments like 
Neutrino Burst Experiment (NuBE) might detect ~20 events per year of ultra high energy 
neutrinos (E >TeV) coinciding with GRB's [ESJ. Among other experiments aimed at studying 
ultra high energy cosmic rays including neutrinos we find the OWL-Airwatch project which 
could detect ~ 3 x 10 3 — 10 5 events (E > 10 20 eV) per year |29 , 30fl . This experiment also 
can look for time correlations between high energy neutrinos and GRB's. Complementary to 
the GRB effect, there is the possibility of looking for neutrino oscillation effect as induced by 
quantum gravity H . An analysis similar to that performed in the case of atmospheric neutrino 



oscillations [§, |31| , can help to set bounds on the parameters entering the description as we will 
suggest below. 

In summary, astrophysical observations of photons, neutrinos or cosmic rays could make 
possible to test quantum gravity effects or at least to restrict quantum gravity theories. 

Indeed, an alternative approach to quantum gravity is based on string theory |3^| . On 
such a basis modified dispersion relations of the type ([I]) have been regained ||33|| . The main 



difference in the case of photons is the absence of helicity dependence that is present in the 
loop quantum gravity approach. As for the case of spin 1/2 particles both approaches seem to 
agree in its helicity independence to leading order of the corresponding effect |Q, ^ |3~5[] . 

Further comments in order here are: delay time effects for traveling particles have also been 
considered on a different basis, for example an open system approach f36|. Also, perturbative 
quantum gravity has been considered to obtain effective dispersion relations |37] . 

Noticeably, the effect considered thusfar involves a Lorentz symmetry violation ||| |4|, |34|, Q , 
which seems to be in agreement with some astrophysical and cosmological scenarios |TT], [13], 
|3~8[ |39| . In this way, these studies naturally overlap with the systematic approach developed 
by Colladay, Kostelecky and collaborators |4(| which provides the most general power count- 
ing renormalizable extension of the standard model that incorporates both Lorentz and CPT 
violations. This framework has been widely used to set experimental bounds upon the inter- 
actions that produce such violations and the observations performed so far cover a wide range 



of experimental settings jlT]. Additional progress in establishing bounds to such symmetry 
violation can be found in Refs. |], [TJ, [4^, [I3[ . 

Finally, it is interesting to emphasize that Planck scale corrections to either particle propaga- 
tion or interactions need not necessarily imply violations of Lorentz covariance jfP| . Recently 
the use of non-linear representations of the Lorentz group, leading to what has been called 
special relativity with two invariant scales, allows for a systematic construction of theories 



exhibiting these features [45 



The paper is organized as follows. In Section 2 we summarize the basics of loop quantum 
gravity in the case of fermions plus gravity. Section 3 is devoted to review the regularization of 
Thiemann for the corresponding Hamiltonian constraint. Section 4 explains general aspects of 
our approximation scheme whereas Section 5 provides the details of the calculations. In Section 
6 the effective Hamiltonian for non- interacting spin 1/2 particles is obtained. In particular, in 
vacuo dispersion relations are given. Section 7 contains some preliminary estimates of the 
parameter T. To conclude, a summary and discussion of our results is presented in Section 8. 
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2 Loop quantum gravity 



This section provides the basic ingredients of this approach, also known as quantum geometry 



| 46| , which we shall use in the sequel. Among the main results along this approach one finds: 
i) well defined geometric operators possessing a discrete spectrum, thus evidencing discreteness 
of space [46fl, ii) a microscopic account for black hole entropy [f|7| and, more recently, hints on 
quantum avoidance of a would be classical cosmological singularity ||8| . (For a review on these 



topics see for example Ref. |49| .) 

To begin with it is assumed that the spacetime manifold M has topology £ x M, with E 
a Riemannian 3-manifold. Here a co-triad e\ is defined, with a,b,c, . . . being spatial tensor 
indices and i,j,k,... being su(2) indices. Thus the corresponding three-metric is given by 
lab = e a^l- In addition, a field K l a is defined by K a b = sgn[det(e{)]K l a e l b) which is related to the 
extrinsic curvature of S . A canonical pair for the gravitational phase space is (K^, Ej/k,), 
where Ef = ^e abc 'eij^e\e k c and k is Newton's constant. It turns out that such a canonical pair 
yields a complicated form for the Hamiltonian constraint of general relativity. A convenient 
canonical pair, making this constraint polynomial, was introduced by Ashtekar []50|1 . Neverthe- 
less, two severe difficulties to proceed with the quantization remained: (i) the implementation 
of a diffeomorphism covariant regularization for the density-weight two Hamiltonian constraint 
hereby obtained and (ii) the extension to non-compact groups of the diffeomorphism covariant 
techniques already developed for gauge theories with compact groups |5Tf| . In fact, the Ashtekar 
variables ( w A l a = T l a — iK l a , iE" 1 / n) |5(J , with T l a being the torsion free connection compatible 
with e*, are complex valued. Namely the gauge group is SL(2,(F), which is non compact. 

Some proposals to come to terms with difficulty (ii) were: to consider real connection 
variables fl52f , to implement a Wick transform [|53| and to define tractable reality constraints 
54]. All of these left open (i). Thiemann subsequently proposed to solve (i) and (ii) by 



incorporating real connection variables while keeping the density weight one character of the 
Hamiltonian constraint. He further provided a quantum version of the theory in the pure 
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gravity well as in those cases including the coupling of matter to gravity ||55|| . His 

approach is next reviewed, since we rely upon it for our analysis of the fermionic case. 

As for the fermionic sector a convenient canonical choice is (£, ir = Here £ = (det g) 1//4 ?7 
is a half density and rj is a Grassman SU (2) spinor. Half densities are convenient because they 
do not lead to cumbersome reality conditions at the quantum level and, furthermore, they 
do not require a complex gravitational connection. However, problems with diffeomorphism 
covariance can emerge. For this reason a further canonical transformation in the fermionic sector 



is necessary in order to dedensitize the Grassman fields through: 9{x) := J s d 3 y^5(x, y)£(y) 
and f (x) := J2 ye s \/K^y)0(y). 

The pieces of the gravity- spin- 1/2 system Hamiltonian constraint read 

#Ei„stein[iV] = / d 3 x N ±= tr (2 [K a , K h ) - F ab ) [E a , E b ]) , 

# S pin-i M = j[ d * x NEi \/M(q) ( mTT{Dai + Va + \ K ^ Ti + CX ) 

+ ^(f^ 2 )Z + * T (^» (3) 

Here f = — |a, where a = {a 1 } are the standard Pauli matrices and we have included an explicit 
mass term. Only one chirality fermion is used 77 = ~(1 + 75)^, which is equivalent to have a 
Dirac spinor ^ T = (r] T , /x T ) satisfying the Majorana condition. The classical configuration space 
is then AjQ of connections modulo gauge transformations, together with that of the fermionic 
field. 



The quantum arena is given as follows [51|. As in any quantum field theory, because of 
the infinite number of degrees of freedom, an enlargement of the classical configuration space 
is required. This is far from trivial since the measures defining the scalar product, which are 
required to provide a Hilbert space, get concentrated on distributional fields and hence lie 
outside the classical configuration space. The key idea to build up such an enlargement is to 
make Wilson loop variables (traces of parallel transport operators) well defined. The resulting 



space AjQ can be thought of as the limit of configuration spaces of lattice gauge theories for 
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all possible floating (i.e. not necessarily rectangular) lattices. Hence, geometric structures on 



lattice configuration space are used to implement a geometric structure on A/Q. This enables 
to define a background independent calculus on it which, in turn, leads to the construction of 
the relevant measures, the Hilbert space and the regulated operators. 

It turns out that the Hilbert space of gauge invariant functions of gravitational and spinor 
fields is given by ]55[ 



ftinv = L 2 ( [A SU (2) ® S\ /g, dfi AL (SU (2)) ® dfi F ) (4) 

Here Q denotes the action of SU(2) on all fields at every point of E. S denotes the infinite 
product measurable space <£) X S X , S x being the Grassmann space at the point x. dfip is the 
corresponding measure on S. An orthonormal basis using \i F can be built up. Let the spinor 
labels be ordered A = 1,2. Take v a finite set of mutually different points. For each v G v denote 
by I v the array (Ai(v) < ■ ■ ■ < Ak(v)), with < k < 2 and Aj(v) = 1,2 for each 1 < j < k. 
Also set \I V \ = k. Thus fermionic vertex functions F^f can be defined as F v f := U. V ^F V>I , 
F v j v := H.j = i9Aj(v)(v ) yielding the desired orthonormal basis. 

Gauge invariant objects are spin network states defined as follows. Take 7 as a piecewise 
analytic graph with edges e and vertices v which are not necessarily closed. Every edge can be 
read as outcoming from a vertex by suitably subdividing edges into halves. Given a connection 
A a we can compute the holonomy along the edge e: h e (A). To each e we also associate a 
spin j e corresponding to an irreducible representation of SU(2). In addition, we attach to each 
vertex an integer label n v , < n v < 2, and a projector p v . For a given n v , one considers the 
vector subspace of Q v (the vector space spanned by holomorphic functions of 9a(v), spanned 
by those vectors Fj >v such that |/| = n(v). The projector p v is a certain SU(2) invariant matrix 
which projects onto one of the linearly independent trivial representations contained in the 
decomposition into irreducible representations of the tensor product consisting of (i) the re- 
fold tensor product of fundamental representations of SU(2) associated with the vector subspace 
of Q v spanned by the Fj jV , |/| = n(v) and (ii) the tensor product of irreducible representations 
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j e where e runs through the subset of edges of 7 starting at v. The resulting gauge invariant 
states are denoted by 

which extend the definition of the matter free case. They form a basis of 7Yi nv - Although not 
orthonormal it can be transformed into one that it is by suitably decomposing the fermionic 
dependence into an orthonormal basis of the Q v . 

To extract physical information we will further need a state describing a flat continuous 
space £ at scales much larger than the Planck length, but not necessarily so at distances 
comparable to Planck length itself. States of this kind were introduced under the name of 
|56f for pure gravity. Flat weave states \W), having a characteristic length C, were first 



weave 



constructed by considering collections of Planck scale circles randomly oriented. If one probes 
distances d » C the continuous flat classical geometry is regained, while for distances d « £ 
the quantum loop structure of space is manifest. In other words, one expects a behavior of 
the type (W^c/aft 1^0 = 5 a b + O (^) . It was soon realized that such states could not yield a 
non trivial volume due to the lack of self intersections ||571| . Couples of circles, intersecting at a 
point, were also considered as specific models of weaves to overcome this defect p8| |. With the 
recent advances on the kinematical Hilbert space H, anx it became clear that all proposed weaves 
were afflicted by two undesirable features. First, they are defined to be peaked at a specific 
(flat or curved) metric, but not at a connection. This is in contrast with standard semiclassical 
states in terms of coherent states, for example. Second, the known weave states do not belong 
either to T~C aux or to a dense subspace of it |53] . It may be possible to come to terms with such 
difficulties by defining coherent states for diffeomorphism covariant gauge theories [^D| or by 
implementing a genuine statistical geometry |JTJ, for instance. Both approaches have recently 
achieved substantial progress. 

Nonetheless, in order to extract physics, there is the alternative possibility of using just the 
main features that semiclassical states should have . Namely, peakedness on both geometry 
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and connection together with the property that they yield well defined expectation values of 
physical operators. An advantage of this alternative is that one may elucidate some physical 
consequences before the full fledged semiclassical analysis is settled down. Indeed, such al- 
ternative may be considered as complementary, in the sense of hinting at possible features of 
semiclassical states which could be further elaborated. After its completion, a rigorous semi- 
classical treatment should tell us whether the results arising from this alternative turn out to 
hold or not. The weakness of the treatment resides on its generality, since no detailed features 
of the would be semiclassical states are used -as opposed, say, to the original weave states- and 
hence a set of numerical coefficients cannot be calculated. Evaluating them will be the task of 
the rigorous semiclassical treatment. 

On top of the purely gravitational semiclassical states, a generalization is required to include 
matter fields. For our analysis it will just suffice to exploit the same aspects of peakedness 
and well defined expectation values, extended to include the case of the fermion field. The 
semiclassical states here considered will describe flat space and a smooth spin-1/2 field living 
in it. Such a state is denoted by \W,£ > and has a characteristic length £. Since no detailed 
information is used on how the semiclassical state is constructed in terms of, say, a graph, 
as opposed to weave states, the present approach yields results relying only on the following 
assumptions: (i) peakedness of the states, (ii) well defined expectation values and iii) existence 
of a coarse-grained expansion involving ratios of the relevant scales of the problem: the Planck 
length £p, the characteristic length £ and the fermion wavelength A. States fulfilling such 
requirements are referred to as would be semiclassical states in the sequel. 



3 The regularization 



We will use the regularization method of Thiemann |x| in the following. We focus on the first 
term of H sp [ n _i/ 2 [N}, Eq. (^). The other terms can be dealt with similarly or, in the case of 
K % a dependence, along the lines of the purely gravitational sector i?Einstein[^] as shown in the 
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sequel. The identity - {A l a , V} = 2sgn(detel)e 1 a yields 

, [N] = * f* xNix) u^im^mm^E^] ^ riDtM 

s P m V 1 2k 2 J K J \ y/q(x) J J 

(6) 

where 5 is small and V(x, 8) = J d 3 y Xs{ x i y)\f q{. x ) is the volume with respect to the metric 
q a b of a box centered at x. xs{x,y) = Y\l=i 9(8/2 — \% a ~ U a \) is the characteristic function of 
the box. Next one changes from £ to 9 variables and introduces 



Y,tf{x){T{D a 6) A (x)6 A \ 



x 



(7) 



which regularizes (|6]) as we show now. Here /° is a real valued &dsu{2) vector field. By defining 
(d a 9 A )(x) '■— lim^o d x a9 A (x), and recalling 



B A (x) = J d s y^5jx~y~)Uy) = ^a 



(8) 



with 9 A = J d 3 y Xc ^^ A, one gets 



(d x aO A )(x) = / d 3 y 



(d y a9 A )(y) 



(9) 



Now let us divide S into boxes of Lebesgue measure e 3 centered at x n and let Eq.(^) be the 
limit e — > of 

^/f(a; n )(r^ £ )A(x n )^(x n ) , (10) 



which upon expressing 6 |e in terms of £ becomes 



d x / cry 



[r i( 9 a £(x) + (^a0rn)£(>))] A 6i • 



Notice that in spite of the density weight of £ no Christoffel connection is needed since it 
will drop after the h.c part is considered. In the limit e — > 0, Xt( x > x n)/t 3 — > S(x,x n ) and 
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Xe(x, x n ) — > Sy tXn , so that after x integration and sum over x n the result is Eq.(||) where ff has 
to be properly identified. 

To proceed with the quantization one keeps fermionic momenta to the right and replaces 9 
by fogf- in (0). By further multiplying and dividing by 5 3 the factor 5 3 ^/g(x) in the denominator 
is changed to V(x, 5). This can be absorbed into the Poisson brackets already appearing in @. 
When Poisson brackets are replaced by commutators times 1/ih one obtains an operator whose 
action on a cylindrical function / 7 associated to a graph 7 and containing fermionic insertions 
9a at the vertices v G V(j) gives 

h 



H 



(i) 

spin— 1/2 



W 7 



2£1 



v£V{~{) x 



X 



X 



A* a (*),V^W) 




V(x,5) 



(T k V c 6) A (v)- 



-$x,v + h.c. 



A- 



(12) 



To complete the regularization an adapted triangulation to 7 of £ is introduced. Using 
h s (0,6)9(s(5) — 9(s(0)) = 5 s a (0)(V a )9(s(0)) and replacing the connection operators by ho- 
lonomy operators allows one to absorb the 5 3 factor. Also one replaces V(v, 5) — > V v . Because 
of the presence of S XiV , becomes concentrated in the vertices ^2 v (A)=v' Hence we obtain 

mp 



H {1) 1/2 

spin— 1/2 



[N] 



2£l 



E E 

v£V(-y) v(A)=v 



AjkJJK 

e J e 



x tr I Tih Sl ( A) 



ra *i(A)> 



tr I Tjh Sj(A) 



x 



[(r k h SKi5) 9)(s K (A)(5)) - 9(v)] 



d 



A d9 A (v] 



+ h.c. 



(13) 



which is the operator we shall use below. 
At this stage we introduce the notation 



Wua 



tr f Tjh 



i'lsj(A) 



s/(A)> 



(14) 
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4 General structure of the calculation 



The effective Dirac Hamiltonian H F is obtained by considering the expectation value of the 
fermionic sector of the quantum Hamiltonian, with respect to \W, £). Inside this expectation 
value, operators are expanded around all relevant vertices of the triangulation in powers of the 
segments s£(A), having lengths of the order lp. In this way, a systematic approximation is 
given involving the scales lp <K C < Xp. Here is De Broglie wavelength of the neutrino. 
Corrections arise at this level. Let us start by taking the would be semiclassical state (WBSC) 
expectation value of Thiemann's regularized Hamiltonian 

h 



E(v) 



e ijk e IJK x 



x{(W,C\ 







de{v) 



d 



-r k 6{v) 



w iIA (v)w jJA (v)\W^) 
Wim(v)w jJA (v)\W,Cl - c.c. } 



(15) 



where 



36{v) 



_d6{v) 

is the fermionic momentum operator . 
Our strategy to use the proposed regularization as a computational tool will be to expand 
the expression ( |15"D around each vertex of the triangulation in powers of the vectors sj(A). To 
proceed with the approximation we think of space as made up of boxes of volume C 3 , whose 
center is denoted by x. Each box contains a large number of vertices of the graph associated 
to the WBSC, but is considered infinitesimal in the scale where the space can be regarded as 
continuous, so that we take C 3 ~ d 3 x. 

Next we discuss how to estimate the average contribution in each box. To begin with we 
consider that the volume of the box is 



d 3 x^ J2 ' 

i>6Box(x) 



Ev) 



(16) 



and define the average < T(x) > of the quantity T(v) defined in each of the vertices contained 
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in the box, as 



d 3 x < T(x) >= 

6 Box (if) 



E v) 



T(v) 



(17) 



The WBSC expectation value of the Hamiltonian (|15|) is of the type 



E(v 



■{W,Z\F(v)G(v)\W,0 



£*•(*) E^ 

Box(x) iiGBox(i?) 



.(w,t\jrG(v)\w,t) 



E(v) 

= [ F{x) < ^ G(x) > d 3 x (18) 

Here, F(v ) is a fermionic operator which produces the slowly varying function F{x) within 
the box; i.e. £ « Xjj. Furthermore, < 4- G(x) > is the box average of the rapidly varying 
WBSC expectation values of the gravitational operator -k-G(v) within the box, whose tensorial 
and Lie-algebra structure is determined by the tensors characterizing the continuum spacetime 
we are dealing with, i.e. e M , r fc , e cde , e fc ' m , where q ab —e m is the corresponding classical 
3-metric. The order of magnitude of these box averaged quantities is estimated according to 
some prescriptions to be specified in the sequel. The method does not provide exact numerical 
coefficients which can only be obtained from a detailed knowledge of the semiclassical state. 

Since we assume the fermionic fields to be slowly-varying functions inside each box, we 
demand the following behavior of the fermionic operators inside the WBSC 



(w,t\....e B (v) 



d 



de A {v) 



\w,o = e 



-4b(v)....--k a (v). 



(19) 



where the normalization constant B is to be determined in such a way that the zeroth order 
approximation reproduces the corresponding classical kinetic energy term in the Hamiltonian. 
In this way we have 



d 



4K 



E E Mi 



. e i i k e IJK x 



d9(v) 



Box(x) d£Box(x) 
Tkh SK (A)0(v + S K (A)) 

d 



_dd{v 



-T k 6{v) 



E(v) 

W iIA (v)WjjA{v)\W,£) 

Wim{v)wjjA{v)\W,Cl - c.c. } 



(20) 
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In the above equation we have written only the fermions in its classical version and we have kept 
them inside the WBSC expectation value just in order to be able to write the above expression 
in a more compact way. The holonomy h SK / A ) still contains the gravitational connection. In 
order to have a convenient bookkeeping of the terms involved, we write Eq.(^UJ) as 



+5 E E^ "<« 

Box(a;) v&Tiox(x) 



Ev) 



_ € ijk e UK x 



a 



p9{v 

x w iIA (v)w jJA (v)\W,£) - c.c.} 



s a K V a (v) )+-( s a K V a (v) ) ( s b K V b (v) ) + .... ) 6(v) 



(21) 



where the derivative V a includes the covariant derivative Va which is compatible with the 
classical metric that we are considering and a piece A l a (v) producing the quantum corrections 
associated with the beginning of the continuum at the scale C . Upon taking WBSC expectation 
value we can make the following substitution 

Va{v)9{v) — fva^+i^hW) (22) 

From now on we restrict to a continuum flat metric, in such a way that 

Va (X) = a . (23) 

The partial derivative does not change inside each box, which is consistent with the idea that 
each box is supposed to represent a point of the continuum. In this way 



iQ 



x {(W^\ 



E E N ( 

Box(a;) fSBox(x) 

d 



E(v) 



(W{i;) r k [s a K V a {v) + ^s a K s h K V a {v)\ h (r) \ 



9(v) 



x WiiA(v)w jJA (v)\W,£) - c.c.} 

Our problem now is to parameterize and to estimate the general spinor 

Y? n) (v) = c y V^(W,e|^ fl S....*V ai (t;)V a9 (t;)...V aB (t;)e(t;) 
xw iIA (v)w jJA (v)\W,£) 



(24) 



(25) 
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in a given box. Here n denotes the number of covariant derivatives in the expression. Using 
the above we can write 

RF = + W £ £ £ Uv)^)n(v)r k Y^v)-c.c. (26) 

P n=l,2,... Box(x) «eBox(f) V V 

Before calculating the corresponding box averaged contributions let us obtain the expressions 
for some of the first quantities YFJv). A direct calculation shows 

Y(i)(v) = s a K (v) (X kK (v)d a m + X? kK (v)r m ) (27) 

Yfo(v) = s a K (v)s b K (v)X kK (v)d a d b ax) + 
+s a K (v)s b K (v)X™ kK (v)r m 2d b ax) 

+s a K (v)s b K (v)XT kK (v)T m r n ax) } (28) 

Y(%(v) = s a K s h K s c K (X kK (v)d a d b d c + 3X™ kK (v)T m d b d c +)^) 
+s a KS b K s C K (3X a T kK (v)r m r n d c + X^ kK (v)r m r n r p ) 

where we have explicitly enforced the symmetry implied by the factors s a K s h K s c K .... in Eq . (|25|) . 
The gravitational quantities, which are rapidly varying inside the box are 

(W,£\w iIA (v)w jJA (v)\W,0 
e«"e iJK (W,t\A™(v)wuA(v)w jJA (v)\W,0 

(W^\A^(v)A^(v)w UA (v)w jJA (v)\W,0 
(W,^(v)A^(v)A^wua(v)w 3 j A (v)\W,0 
Next we write the corresponding general expressions 

l=n,n-l,....l,0 

x r mi+1 r mi+2 . ..r mn d ai d ai ....d ai £(x) (29) 
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X kK \ 




AjkJJK 






AjkJJK 

= e J e 


\rmn kK , 
A ab 1 




AjkJJK 


yinnj) kK 
^abc 1 




AjkJJK 

= e J e 



with 

KlXZtr.2 nkK (v) = <*V J *<w;ei4£i(T;)^ 

(30) 

In the above expressions the partition (1,2,.../) and (I + 1, I + 2, ....n ) is made in such 
a way that the indices of the first set count the number of partial derivatives acting upon 
the fermion. The indices of the second set count the number of r-matrices involved which is 
equal to the number of additional connections, i.e. besides those contained in the combination 
wuA(v)wjj&(v), appearing in Eq.(|30|). 

The box-averaged expression for the Hamiltonian reduces then to 

RF = Z§- E fd 3 x7r(x)r k <Y^(x)>-c.c. (31) 

n =l,2 J 



In more detail this is 



rF = + wr E / rf3f n W Tk x 



P n=l,2,.. 



2=n,n-l,....l,0 

r m i+1 r m; + 2 • ••' r m„ ^aj ^ai • • • ■ £ (#) ■ (32) 

Here it is convenient to define the auxiliary quantities 

rpa 1 a2....aimi +1 m l+ 2 m n k _ ® ^ / Oj Oa „ a ; \ / a i+i a i+2 o„\ j^m i+1 m( +2 ..m„ fc_ftT / \ ^ 

^£7 V X A""" A"/ i K K K ) a i+1 a ;+2 ....a n V / 

(33) 

in such a way that 

RF = E E (") [d 3 xin(x)r k r*<»"" a ™+™+* mn r mi+1 T mi+2 ...r mn 

xd ai d ai ....d ai £(x) 

(34) 
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n=l,2,... l=n,n— 



RF =- E E (35) 

n=l,2,... i=n,n-l,....l,0 

The box-averages < F(v ) > are subsequently estimated using the corresponding dimensions in 
terms of the available tensors in flat space 

°*a = ( ^ r fe ; ^ £C£ fe gfcbn ? afc = gab ^ gij ^ $a (gg) 

in a manner described in the next section. In this way we are demanding covariance under 
rotations at the scale C. As the final input in our prescription we impose that when aver- 
aging inside each box, the order of magnitude of the corresponding expectation values of the 
gravitational operators are estimated according to 

(W,£\...,A ia , ...,VK,---|W,0~--4 (^) <■■■ 'V (37) 

respectively. 

In previous work [[| we have set T = on the basis that the coarse graining approximation 
does not allow for the continuum connection to be probed below ^. On the other hand, by 
adopting say kinematical coherent states for representing semiclassical states, one would set 
T = 1 for two reasons: first, in the limit h — > (£p — > 0) ( p7|) yields just zero, in agreement 
with a flat connection, and, second, such an scaling would saturate the Heisenberg uncertainty 
relation. Nonetheless, physical semiclassical states may lead to a leading order contribution 
with T / 0,1, thus we leave it open here. Now we have all the ingredients to perform the 
expansion in powers of s. 

Let us consider now the contribution arising from the gravitational operators contained in 
wu/\{v) of Eq. (p0|) . We choose the parameterization 

Wha(v) = s/ a w ia + s 7 a s/ w iab + ... (38) 

leading to 

wua(v) Wjj A (v) = s/ a sj b w ia w jb + sj b sf w ia w jbc + sf sf sf w iab w jc + 0(s 4 ), (39) 

18 



where Wi a and Wi a b, which are independent of s, need to be calculated explicitly. The product 
of w starts quadratically in s. 

To this end we will need the r-algebra, recalling that = — with a k being the standard 
Pauli matrices. We have 

1 1 

1 _1 1 1 

ab ~ 2 eabmTm ~ Tc 

TkT a T k = ~T a (40) 

From the definition fli"4]) we obtain 

wua = Tr (ri(VV - h B1 VVhjty (41) 
Up to second order in s, we have 



w 



ja = -Tr(r i ([A I ,VV} + ^[A I ,[A I ,VV}} + ...)] , A T = s T a A ta r l 

= Sl a ~ [A ia , W] + sj a sj h \e ikl [A ka , [A lb , W}] + . . . , (42) 

Z o 



which finally allows to identify 



w 



= -[A^, ^V], w iab = ~e m [A^ [Aw, W]]. (43) 

Z o 

The scaling properties under the semiclassical expectation value of the above gravitational 
operators is 

p 3/2 / p \ ™ T 

(Wt\...^... an ...\Wt)^^U) . (44) 
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5 The calculation 



The detailed calculation of the correction terms is performed according to the following prescrip- 
tion: first we set T = 0, then we consider an expansion to order i 2 p and finally we reintroduce 
the non-zero value for T in the corresponding terms. 
Let us recall the general expressions from last section 

r r a 1 a 2 ....a l m l+1 m l+2 m n k _ ® ^ I a x „a 2 »[\ f"l+l„ n l+2 a n \ 

1 ~ ^£7 ^ \ b K b K--- b K) \ b K K ■■■ b K) 

x x mi+ r i+2 'r n kK (v) > (46) 

xA^(v)w iIA (v)w jJA (v)\W,0 (47) 

The partition is made of the number n such that we have / derivatives and (n — I) r-matrices, 
all arising from the presence of n covariant derivatives in Eq. (|24"D . We recall that 



Wua{v) w jJA (v) = (s 7 a s/ w ia w jb ) + (s/ a s/ s/ w ia w jbc + s/ a Sj 6 sj c w iab w jc ) 

+0(s 4 ), 

= (wua(v) w jJA (v)) 2 + (w iIA (v) w jJA (v)) 3 + 

- + {wua(v) w jJA (v)) N + (48) 

where N counts the powers of s in the term wu A {v) Wjj A {v). Under the semiclassical state we 
can estimate the contribution of each term in Eq. (f48|) as 



{W,£\ (w i iA(v)w jJA (v)) N ....\W,0=£ 3 P ( ^ ) , (49) 



e_p\ 



7V(1+T) 



where we have used Eqs. (|37D and (f44|). 
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5.1 The case n = 1 

Here we have two possibilities which we consider separately: (n — 1, I — 1) and (n — 1, I — 0). 



5.1.1 Case n = 1, / = 1 

K( 1= ij d 3 xn(x) r k T« fc c> a £(£) (50) 

= (51) 

= e^V J ^(H/,e|^/A(^)^jA(^)|^,0 (52) 

We present this first calculation in some detail since it sets the stage for all the remaining 
estimates 

T"* = ^ <s a K e^ k e IJK (W^\w UA (v)w jJA (v)\W,0 > 



ft ff 

E «™4 (|) 



ft / ft \2( 1 + T ) //> \3(1+T) 



^ = ^(|) 2(1+T, (, m + , 113 (|) <1+r, ....| (53) 

The notation is i? n jjv and these numbers are assumed to be of order one. Choosing i? 112 = 1 we 
are left with 

ft ff \ 2 (!+ T ) 

T " = s °*w P \i) (54) 

where 



Fn= I l + ^iia l^-J ....)• (55) 
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Our convention in the sequel is to write T n i so that its first term is $ n ;2, which is a pure number 
not depending either on or on £. In this way 

H n = TFT HH ^ii / d3 % i7[ ^) r a d a £(x) (56) 



In order to recover the standard kinetic term we have to choose 



r \ 2(1+T) 

9 = 44 ( - , (57) 



which reduces to our choice in Ref. [|J, for T = 0. We obtain 

Hfi = Fn y d 3 x iV(f) m{x)T a d a i{x) 

as the final result for the kinetic term piece of the Hamiltonian 

5.1.2 Casera=l,Z = 
The basic quantities are 



Under the scaling properties we obtain 



T rnk = S mk}_ll_P} ^ (G;|) 



T 



with 

(l+T) 



Hf = -\f*> J d 3 fz7r(£)£(f) 



(58) 



Hf = 1 J d 3 x tt(£) r k r m T™ fc £(£) (59) 
T i" fc = Tin < s k e ijk e IJK (W^\A^(v)AiA(v)w jJA (v)\W^) > (60) 



Fio = (^102 +^103 (^) +-.. ) (62) 



(63) 
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5.2 The case n = 2 



Here we have three possibilities: (n — 2, 1 — 0), (n = 2, 1 — 1) and (n = 2, / = 2) 

5.2.1 Case (n = 2, 1 = 2) 

H& = i y rf 3 x 7r(f) r fc T 2 f fe 9 a ^(f) (64) 

* =^T< s k s k e l3 h IJK (W^\w tI A(v)w jJA (v)\W,0 > (65) 

From the above equation we see that T 2 a | k must be symmetric in the indices a, b. Nevertheless, 
the only three index tensor at our disposal in flat space is e abk so we must conclude that this 
contribution is zero, that is to say H^ = 0. 

5.2.2 Case (n = 2, 1 = 1) 

H£ = 2i J d 3 x tt(£) T k T m T£™ k d a Z(x) (66) 

T 2 T " = ^ < ^4 e^V^(W,e|irW^/AW^vA(^)|W,0 > . (67) 
Here we have no symmetry requirement, so the antisymmetric tensor is allowed and we have 

T 2 T h = e amk [^Y F 21 , JF 21 = (tf 212 + $213 + +-.. 1 (68) 

The final contribution is 

l+T 



H 21 



-2JT 21 J Sxm(x)r a d a i(x) (^j (69) 



This is a correction of order (^) to the standard kinetic term. 
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5.2.3 Case n = 2, I = 

Hf = i y rf 3 f 7r(f) r fc T 2 " m fc r n r m e(f ) (70) 

^"o mfe =§r< s kSk t l3k t IJK {W^\A-{v)A™{v)w UA {v)w 3jA {v)\W,i) (v) > (71) 

Here we do have symmetry requirements. Since the operators A™{y) and A™(v) commute, the 
above tensor must be symmetrical in the indices n, m. Again, in flat space the only three index 
tensor at our disposal is the e nmk , which leads to Hf = 0. 

5.3 The case n = 3 

Here we have four possibilities: (n — 3, I — 0), (n — 3, I — 1), (n — 3, 1 — 2) and (ra = 3, / = 
3). 

5.3.1 Case (n = 3, Z = 3) 

Hf 3 = i y rf 3 f ir(x) r k T 3 f k d a d b d£{x) (72) 

T 3 f k =^r< ^4^^V^(w,e|^AW^vA^)|^0 > (73) 

Here we have to impose the symmetry in the indices a, b, c. In flat space the required tensor is 

tf 3 ck = - (S ab 5 ck + 5 bc 5 ak + 5 ca 5 bk ) (74) 
3 

which turns out to be symmetric in all four indices Then 

iabck ±abck fi2 <r- t I „Q i „Q / " 



T 3 f = ^33, ^33 = | ^332 + ^333 ( f ) + - ) (75) 



Combining the above results we obtain 

H 3 F 3 = F 33 Jd 3 xtn(x)T a d a d 2 ax)£ 2 P , 
the leading term of which is independent of the scale C. 
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(76) 



5.3.2 Case n = 3, 1 = 2 

Hf 2 = 3y ftm(i) r fc T 3 f™ V m d a c> b £(£) (77) 

^ mfe = -gr < ****** e ijk e IJK (W,^(v)w a A(v)w jJA (v)\W,0 > (78) 
The above tensor must be symmetrical in the indices a, b only. In flat space we can construct 

tf 2 mk = a 32 S ab 6 mk + f3 32 (5 am 5 bk + S bm 5 ak ) (79) 

Then we have 

1 // \ T / // \ (1+Y) \ 

T abm k = t <gnkf p _ ( _Z \ ^ ^ = I ^ + ^ (£. \ + .... 1 (80 ) 

which leads to 

Hf 2 = "3 (j«32 + ?32 J d 3 X iTT(x) fpd^x) I [^j (81) 

5.3.3 Case n = 3, 1 = 1 

H£ = 3 J d 3 xin(x) r k T*r k r m r n d a m (82) 

T 3 T nk = < ****** c°V JX (w;$|ir(«)^(«)^A(t;)T&iJA(t;)|w;o > (83) 

The above tensor must be symmetrical in m, n. We have 

t arank = ^gmngak + ^ (gmagnk + finafirnkj ^ 

Then 

ff \2(1+T) / x(l+T) \ 

T arank^ t arnnk\P\ ^ ^31 = Ul2 + 1>313 (f J + .... J (85) 

which produces 

Hfi = -^3i (9a 3 i + 6/? 3 i) y d 3 a? i7r(x)r o a a e(f) ff J (86) 
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5.3.4 Case n = 3, / = 

Hf = J d 3 x iit{x)T™ np k T k T m r n r p m (87) 

TT pk =^r< ^4^V J ^,£|i^)^ > (88) 

The above tensor is symmetrical with respect to the indices m,n,p. This means that we need 

t mnpk = 1 (gmngpk + 5 np § mk + gpmgnkj (gg) 

so that 

// \3T / , f x(l+T) \ 
T 3 7 P ^ = tZ^f, ffl ^30, ^30 = Uo2 + ^303 (j) + J , (90) 

which produces 

Hf = ^30 / <?x iv{x) m f (I) 3T (91) 

5.4 The full contribution H F 

It is given by 

H F = £ (92) 

where the different pieces were calculated in the previous subsection. Next we make a clear 
separation among terms containing either £ P or C and those purely numerical factors not 
including these quantities, through the following series of redefinitions. First we change T n \ 
into Q n i, where 

/ f£p\ {1+r] \ 

Qnl = ( Xnl2 + Xnl3 ( J + ) , (93) 

by absorbing all numerical factors in the latter coefficients, in such a way that the Hamiltonian 
is written as 

r ( ft \ 1+T f p \ 2 ( 1+T ) 

H F = J d 3 x N(x) m{x) Gn + G 21 | + &i ( ^ J + fe^V 2 + 
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xr a <9 a £(f) 

+jsx i tt(x) i f f? 10 + + ^4v 2 + ..J e(x) 

In particular we have Qw = and what we have done amounts effectively to a numerical 
redefinition of each finiN into >i n iN- Finally we factor out the contributions arising from different 
powers of (^) and ip. To the order considered we have 

fO \ 1+T tP \2( 1 + T ) tP \ 1+T \2( 1 + T ) 

6u + &1 (|) +631 (|) + ... . 1 + K1 (|) +K2 (|) + ... 

= f (94) 

where Kp are numbers not containing either ip or C. The factor one in the RHS of the first 
equation arises from the condition that in the limit ip — > we recover the standard fermionic 
Hamiltonian. Analogously we redefine 

/| \2(1+T) \1+T /| \2(1+Tj 

£l0 + ^30 f J = K-4 + K-5 f ~£ J +Ke [^) + 

G32 = y (95) 
where the remaining k's fulfill the same condition as the previous ones. The above leads to 

H F = J d s xN(x) m{x) ^1 + Kl (y) 1 ^ + ^ (y) 2(1+T) + y^V 2 + - ) r a d a m 

1 /fi \ r / /I) \ l+T //? \ 2(1+T) 



(96) 

which is our final expression for this piece of the effective Hamiltonian. 
5.5 The mass term 

Since the spinors £, n are half-densities regularization of the mass term requires a treatment 
along the lines of |55[ that effectively dedensitize them. For our purposes, however, it is enough 
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to determine the leading contribution and this can be done consistently with our strategy above. 
Namely, using the triangulation again, 

= ^AT(^(A))^^6^6^(A)4(A)^(A)e T ( S (A)) (^ 2 ) h s{A) Z(v(A)) + c.c. 

(97) 

and in the quantum theory, by adapting the triangulation to the flat WBSC state with fermions, 
we get 

E m = (W,Z\H m \W,0 = ^AT( v (A))^|^e"V bc 4(A)4(A)^(A) 

x(W^\9 T (s(A)) (ia 2 ) h s(A) 9(v(A))\W,0 + cc 
= Yl ^(^(A))^ y 6^6^(A)4(A)^(A) 

x((W,Z\6 T (v)(ta 2 )6(v)\W,0 

+ (W,Z\s a d a 6 T (v) (ia 2 )6(v)\W,0 
+(W^\9 T (v) (ia 2 ) s a A a (v) 6(v)\W,0 

+ (W,£\s a d a 9 T (v) (ia 2 ) s a A a (v) 9(v)\W,0 + -- - ) +c.c. 

(98) 

Hence, to leading order, the modifications coming from Planck scale to standard flat space mass 
term for fermions are just 



H m = 'A(f( W 2 )( + K9 f P f (ia 2 )r a d a ^ 



( |) e^ 2 )i + «ii ip e («r 2 ) r a d a £ + c.c. ] , (99) 



where we have set N = 1. 
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5.6 The extra contributions 

Here we study the second term of Eq. (|3|) 

, : = [ d 3 x^[V a (n T T^)+c.c.} 



(100) 



Repeating the procedure applied to H^^_i step by step with the pertinent modification in the 
fermionic term leads to 



spin- 1/2 L 



JjkJJK 



wGV(7) v(A)=v 

x tr | Tih Sl ( A ) 



h' 1 



tr I Tjh Sj(A) 



X 



(h;i {A) T k h SKiA) 9)(s K (A)(5)) 
d 



d 



a d6 A (s K (A)(5)) 



JA de A (vj 

The corresponding expectation value behaves as 



+ h.c. 



H (2) :) 

spin— i' 



cfa;[Z 1 <9 a £ a + Z 2 ,9 a <9 fe E ab + ... 



(101) 



(102) 



where the Zj's represent the estimate of the non fermionic factors in ( |101| ) and the B's represent 
the fermionic bilinears. Clearly every B term, being a boundary term, will not affect the effective 
dynamics of the spin-^ particle. 

We proceed to describe the third term in Eq. (§) 

r(3) 



spin- 



2^2' 



(103) 



To regularize it we can treat the ^= factor similarly as with H^^_ x term. As for the K % a factor 



we recall 151) 



H E 



-~{Ai,{V,H E }} 

- [ d 3 xe abc tr(F ab {A c ,V}) . 



(104) 
(105) 
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Upon regularization of the quantum H E one gets 



Hi 



E 

A ) 



J^T^r (h atHA) h Sk(A) 



h^ 1 V 



(106) 
(107) 



Following the standard procedure to incorporate holonomies and fermions one arrives at the 
quantum version of (|103|) as 



h 



JJK C LMN C ... ... 



v(A)=v v(A')=v 



x tr ( r fc /i SK lK,tr ( h n h 



xU«(A))- 



5 



+ h.c. \ / 7 



' sjv(A')' Vl> 



(106 



0fc(«(A)) 

where the prefactor was obtained from (^2) (/i) (-^-) ^2^2 1 with the k factors coming 
from the different Poisson brackets identities, the H factors from the quantization procedure 
and the £ 2 P factor arises similarly as in previous steps. We want to estimate the leading order 
contribution of ( |103| ) to the effective dynamics of the spin-| particle. For simplicity in the 
discussion we set T = 0, but an extension can be given for T > 0. The leading order is 
estimated by noticing that under < (...) > the contribution of an holonomy h is 1 if it is 
not an entry of a commutator, otherwise it scales like £p/C Such a leading order contribution 
behaves then like (fJ, 2l -j^J (^f ''£% = £ (^) 4 = m P (^) 4 , which is highly suppressed 
given the quartic order and tp <K C 

6 Effective dynamics for spin-| particles 

The total effective Hamiltonian for our spin-^ particle is then 



H T = H F + H m 



d 3 x 



i 7i(x)T d d d A + c.c. 
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Til 

+ ^(x) {ia 2 ){a + 2hPr a d a )m 



m 



+ ^vr T (f) (a + 2hf3 r a d a ) (ia 2 )Ti(x) 
2h 



(109) 



where 



A | 1 + k x ( ^ 



T+l 



2T+2 



«a ( ^ ) + ? 4 V 2 



C = n I «4 ( ^ 



n = | 1 + Kg ( — 



P 

c 

2T+1 



3T+2 / \ X 

+ f (z) 



T+l 



;no) 



It can be seen that the terms proportional to K4, K5 and kq correspond to a renormalization 
of the fermion mass. They will not be considered in the sequel, although they would give rise 
to standard neutrino oscillations even if the neutrino bare mass is zero. 

The wave equation becomes 



■^9 . ^ * „ C 

ih— - 1 h A a • V H 

dt 2C 

d C 

ih— +i h A a ■ V 

dt 2L 



f (t, x)+m(a- (3 ih a • V) x{t, %) = 0, (111) 
X(t, x) +m(a- (5 ih o ■ V) f (t, x) = 0, (112) 



with x(t, x) = i cr 2 £* (£,#). Following the standard steps one can verify the consistency of the 
above equations. 

Eliminating x from QUI]) we obtain 



1 



X = — 



m a — (3 i ha • V 



d C 
ih—-i hAa-V + — 
at 2C 



(113) 
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Substituting in (|112[ ) we obtain the second order equation 



d C 
ih— + i hAa-V - — 
at 2C 



d C 
ih— -thAa-V + — 
at 2C 



m 2 (a- (3 i ha- V) 2 £(t, f) = 0. 



(114) 



The above equation has positive and negative energy solutions 



W(p, h)e-i Et+ ^ g , W(p, h)e^ Et --^ 



lis) 



where it is convenient to take the spinorial part W(p, h) as helicity (a ■ p) eigenstates, with 
h — ±1, so that 



cos 



The dispersion relation has the following form 



-e ^ sin( 
cos(|) 



;ii6) 



E± = \/U 2 + r 2 )\,7\ 2 ■■: ,n 2 n 2 ( — ) ± C |/7j. 



(117) 



with 



A = | 1 + Kl ( ^) T+1 + k 2 (^) 2T+2 " y 4 I , ^ = ^ ( ^ + 2a(3m' 



C 



2 

a = I 1 + k 8 



T 



T+l\ / „ \ T+l 

P1 ^ ?=&+Xt) ■ 



'1181 



The ± in Eq. ( |117| ) refers to the dispersion relation of the helicities ± respectively. Let us 
emphasize that the solution £(£,x) to either Eq. ( |111|) or Eq. (|112|) is given by an appropriate 
linear combination of plane waves and helicity eigenstates. This is not unexpected since we are 
dealing with massive neutrinos. 
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In the sequel we write down the dispersion relation in terms of an expansion up to second 
order in (^-) T - The coefficient of each power is subsequently expanded in powers of m in the 
combinations either or (m£ P ). In this way we obtain 



E±(p,£) = 



p + ^-±£ P Q™ 2 «gj + 4 (~^3P 3 + ^ (2k 3 + Kg) 



1 £ 



T+l 



'•:/' I ± ; /' 



2 2 

p m 



+ 



2T+2 



where the previous coefficients denoted by k appear in the following combinations 
©ii = (2/ci - Ak 8 ) , ©i2 = (8«n + 2k 7 + 8k 8 k 9 ) , 



(119) 



(120) 
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— 32^1 + 32k 2 + 64kik 8 



(121) 



The velocity is defined as 



v±(p,C) = 



dE ± (p,C) 
dp 



and is 



v±{p,C) = 



m 
V 

T+l r 



+ 4 



~K 3 p 2 + - (2k 3 + rejj) 



m 



+ 



«1 + 



2T+2 



@ n m 2 

o 

m „ 
«2 + — 



=F %(t PP ) 



Up 2 



(122) 



(123) 



within the same approximation. 

We will be mainly interested in the case of ultrarelativistic cosmological neutrinos, which 
mass we take to be m = 10 _9 GeV, in the range of momenta 10 5 < p < 10 10 GeV, where 



p » m, 



m 



{pip) » (me P )-, (pip) 2 » (m£ P y 
P 



(124) 
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Such regime allows us to introduce some simplifications in the coefficients of (^f) T - The results 



are 



m . (\ ,„ , 2 \ f 3 



£p£ ± (p, £) = (p£ P ) + (f P m) — ± - (£ P mY k 9 J - -« 3 (£ P p) 



+ 



T+l 



«1 (Kp) =F «7 



+ 



c 



2T+2 



(125) 



together with 



v±(p,C) = 1 - — - -k 3 (« P p) + 



2p 2 2"°'"^ ' V C 



T+l 



«i =F y (^pP) 



+ 



2T+2 



«2. (126) 



7 Parameter estimates 



7.1 The scale £ 



In order to estimate some numerical values for the modifications to the velocity of propagation 
we must further fix the value of the scale C. Two distinguished cases arise 

7.1.1 The mobile scale £ = A = - 



Recall that £ is a macroscopically length scale, being defined by the given geometry which 
indicates where the non perturbative states of the spin-network can be approximated by the 
classical flat metric. The neutrino is characterized by energies which probe to distances of 
order A. To be described by a classical continuous flat geometry equation, as the found Dirac 
equation, it is necessary that one follows in the correct range with respect to C, i.e. C < A. We 
take the marginal choice of the equality in order to be able to make some further estimates. 
Previous expressions reduce to 



£ P E ± (p,C)\ c __ 



m ■ 1 <- -2 1 ,„ x3 



(p£ P ) + (£ P m) — ± - (£ P m) Kg- -k 3 (£ P p) 



+ W 



T+l 



«i iptp) =F «7 



(ippf 



+ (p£ P ) 2T+2 K 2 (pip) 
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(127) 



and 



v±(p, C)\ c=lp = 1 - ^ - |k 3 {i P pf + (£ pP ) r+1 /mT 1_(7 /(/ ,) +(i r ,,f< + *K, (128) 



2p 2 2 



2 



In the case T = we recover the expressions (9) and (10) of Ref. 
7.1.2 A universal scale for C 



Recently, in [|13| a universal value for C was considered in the context of the GZK anomaly. The 
study of the different reactions involved produces a preferred bound on C : 4.6 x 10~ 8 GeV~ 1 > 
£ > 8.3 x 10~ 9 GeV A_1 . Actually, in this case the general expressions ( |125| ) and ( |126j ) are valid. 

7.2 The exponent T 

In this section we follow the conclusion derived fom the Super-Kamiokande atmospheric neu- 
trino experiment that neutrino oscillations are well described by the mass differences contri- 
bution to the energy (n = —0.9 ± 0.4 at 90% CL, in the standard notation |62]| ) ||] .That is 



to say, any additional contribution to the oscillation length must be highly suppressed. This 
condition will set a lower bound to T which we will use in further estimates. 
The oscillation length L = j^j^ is 

111 1 4nE 

+ W ~ l + fe ' ~ (Amp (129) 

where we have used p ~ E, tp = 1/Mp and Lqg is calculated according to each specific 
additional term in |AJ5|. Since the dominant contribution to the oscillation arises form the 
mass term, i.e. L w L m we must have the condition 



J QG 



>L m ^X, (130) 
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where X is the distance travelled by the neutrinos in each experiment. Let us observe that the 
energy range for current neutrino observatories lies between 10~ 2 GeV and 10 +2 GeV, which 
amply fulfills the condition (p£) < 1. In particular we will consider the case of the SNO 
experiment characterized by the following parameters 



E = 10~ 2 GeV, X = 10 s km = 10 27 l/GeV (131) 

The estimates proceed by considering first the separate contribution of each arbitrary parameter 
Ak and subsequently the case where we have more than one non-zero contribution to the 
oscillation. 

7.2.1 The case £ = A = - 

p 

(i) Ak 3 ^ 0, A«j = Ak 7 = 

The quantum gravity contribution to the oscillation length is 

4tt fM\ 2 1 



l «° = a^{e) e' (132) 

from where we get that 

Ak 3 < 6.28 x 10 17 , 

which is a very weak bound on the variable. The reason for this is that this term corresponds 
to a second order correction in lp. 

(ii) A«i ^ 0, Ak 3 = Ak 7 = 0. 

Here we have 

2tt /M\ T+1 1 



Using A«i ~ 1, we obtain 

T > 0.152. (134) 

(iii) Ak 7 ^ 0, A«! = Ak 3 = 0. 
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Here we have 



87T 



M\ T+2 1 



(135) 



Ak 7 \ E J E 
Using Ak 7 ~ 1, we obtain 

T > -0.87. (136) 

(iv) When two of the Ani do not simultaneously vanish, we get the following situations: 
(iv-1) Aki ^ 

In this case the term proportional to Ak 7 is suppressed with respect to the one proportional 
to Ak\. The comparison among the terms proportional to Ak\ and Ak 3 leads to the threshold 
T = 1. When T > 1 the term proportional to Ak 3 dominates, while the term proportional to 
Aki dominates in the other situation. 

(iv-2) A ki = 

In this case the competition is among the terms proportional to Ak 3 and Ak 7 . The threshold 
here is T = 0, so that we consider only the case T > 0, where the term proportional to Ak 3 
dominates. 

7.2.2 A universal scale for C 

In this estimate we will assume that L = 10 _8 Gef _1 [E3|. 



(i) Ak 3 ^ 0, A«i = Ak 7 = 

Since this contribution does not depend upon C, we obtain the same result as the corre- 
sponding case in the mobile scale. 

(ii) A«! ^ 0, Ak 3 = Ak 7 = 0. 
Here we have 



L « G =A^UJ E- (137) 



Using Aki ~ 1, we obtain 

(hi) Ak 7 7^ 0, A«i = Ak 3 = 0. 



2tt I 
Vp) e 

T > 1.2. (138) 
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Here we have 



?tt ( C \ T+1 M P 1 



Ak 7 \£ p ) EE 
Assuming Ak 7 ~ 1, we get a bound on T: 



(139) 



T > -0.764, (140) 

which is also a weak bound on the variable. 

(iv) When two of the A/«j do not simultaneously vanish, we get the following situations: 
(iv-1) Ak\ ^ 0, which leads to the threshold T > -jf. When T > 17/11 the term Ak% 
dominates and we are back to the case (i). On the other hand, when Y < ^ we obtain the 
interval > T > 1.2. It should be stressed that the latter bound on T is also compatible 
with the bound on the maximum speed of Ref . , which in our case reads 



v±(p,C) ~ 1 - 7TT 



2 

m 



M(f ) T+1 <KT 22 . (141) 



V, 

(iv-2) Are a = 

The threshold here is T = 3/11. When T > 3/11 the term proportional to Ak^ dominates 
over the one proportional to Ak 7 

So, the present data on neutrino oscillations do not proscribe the theory considered in this 
work with a universal scale £ ~ 10~ 8 GeV~ x 

8 Summary and Discussion 

In this work we have derived an effective Hamiltonian exhibiting Planck scale corrections with 
respect to standard propagation for the theory describing spin 1/2 fermions, using an heuristical 
approach based upon Thiemann's regularization within the framework of loop quantum gravity. 
Corrections arise as a consequence of the discrete nature of space which are manifest at Planck 



scale. The effective spin 1/2 particle Hamiltonian given in Eqs.(10£,llC) was obtained by 



taking the expectation value of the regularized version of the quantum operator corresponding 
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to Eq.(|3|) with respect to a would be semiclassical state \ W, £) describing a large scale flat metric 
together with a slowly varying classical spinor field. Only the basic properties of such state 
were used: (i) peakedness both in a flat space metric together with a flat connection for large 
distances d » C » £p, where C can be thought of as the scale that settles the transition 
between a discrete and a continuous description of space, (ii) well defined expectation values, 
(iii) existence of a coarse-grained expansion involving ratios of the relevant scales of the problem: 
the Planck length £p, the characteristic length £ of the state and the de Broglie wavelength 
A of the spin 1/2 particle and (iv) invariance under rotations at scales larger than C, which 
amounts to express the box-averaged values introduced after Eq.fll8|) only in terms of flat space 
tensors. 

The effective theory violates Lorentz invariance and, in analogy with the photon case ||, 
we assume that the effective Hamiltonian so found corresponds to that in the particular frame 
of reference where the cosmic radiation background looks isotropic. 

Some improvements with respect to our original presentation Q are the following: (i) in 
section 5.6 we have elucidated the contribution of the derivative term and the extrinsic curvature 
dependent term in Eq.(3) by showing that they are highly suppressed in powers of £p. (ii) we 
have extended the corrections to the scaling of the connection by including the new parameter 
T, already considered in our discussion of photons ||, in the form 

(W,£\A ia \W, = + 2^) > ( 142 ) 

where T > can be any real number. 

In Section 8 we have estimated some bounds for T, based on the observation that atmo- 
spheric neutrino oscillations at average energies of the order 10~ 2 — 10 2 GeV are dominated by 



the corresponding mass differences via the oscillation length L m in Eq. ( |129|) . This means that 
additional contributions to the oscillation length, in particular the quantum gravity correction 
Lqg, should satisfy Lqg > L m , which is used to set a lower bound upon T. Within the proposed 
two different ways of estimating the scale C of the process we obtain: (i) T > 0.15 when L is 
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considered as a mobile scale and it is estimated by 1/E and (ii) 1.2 < T when the scale L takes 
the universal value C = 10~ 8 l/GeV, according to Ref. |Tj|. In the above estimates we have 
satisfied the condition (pC) < 1. 

Let us observe that according to Eq.( p.27| ) the mass-difference contribution to the neutrino 
oscillation process will be highly suppressed at high energies while other mechanisms, like 
those arising in the quantum gravity framework, could be the dominant ones Ifnfl . It is an 
experimental issue to settle this question. Since the bound upon T strongly depends on the 
dominant mechanism there is the possibility that the exponent T be energy-dependent. This 
situation is not considered in the present approach. Nevertheless, in order to make some 
numerical estimates related to cosmological neutrinos we will make the extrapolation of our 
bounds in T, which have been obtained in the range of a few GeV, to energies of the order 
of 10 5 GeV , corresponding to typical neutrinos arising from Gamma Ray Bursts. Let us 
consider specifically ultrarelativistic neutrinos of mass m = 1CT 9 GeV, energy p pa E = 10 5 
GeV, traveling a cosmological distance L = 10 10 l.y. = 0.5 x 10 42 1/GeV, where pC < 1. Next 
we give an estimate of two types of delay times (for a more realistic calculation including the 
effect of the expanding universe see Refs.|J, |T3|, J3TJ] ) : 

(i) At v pa -y (v — Vo) , which measures the time delay of the arrival of the neutrino coming 
from a distance L, flying with velocity v = with respect to the time of flight corresponding 
to a particle of mass m with velocity vq = -, p ~ 1. 

(ii) At± pa (t>_ — v + ), which is a measure of the birrefringence effects and it is defined by 
the difference in the arrival time of two neutrinos with opposite polarizations 

In the range 0.6 < T < 2.0 and within the two scenarios presented here for dealing with 
the scale C, we obtain the estimates 

10" 10 s < At v < 10~ 4 s, (143) 
At± < 10~ 17 s. (144) 

To conclude we point out that an effective dynamics and dispersion relations for gravity 
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plus matter were thoroughly studied recently by Thiemann and Sahlmann ||21|| , where the 
semiclassical states were taken as coherent states. They included photons and scalar particles 
and their results for the dispersion relations essentially agree with what we obtain here for 
fermions, even the possibility of having non integral powers in £p for the correcting terms, 
which is encoded in our parameter T. Besides they provide a detailed classification of the 
quantum geometry aspects required in defining the semiclassical regime yielding the effective 
dynamics. Yet another avenue to tackle the problem of defining semiclassical states in quantum 
geometry is currently under investigation that establishes a relation between Fock space and the 
kinematical Hilbert space for diffeomorphism covariant theories of connections such as quantum 
geometry/loop quantum gravity [^5], Q (see also [0). It will be interesting to compare all of 
the above proposed semiclassical states in the context of the quantum gravity phenomenology. 
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